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Abstract. We compute the scalar gravitational radiation from a binary pulsar sys- 
tem in the simplest model that exhibits the Vainshtein mechanism. The mechanism 
is successful in screening the effect from scalar fields conformally coupled to matter, 
although gravitational radiation is less suppressed relative to its general relativity pre- 
dictions than static fifth forces effects within the pulsar system. This is due to a 
. ^ \ combination of two effects: firstly the existence of monopole and dipole radiation; sec- 

ondly the Vainshtein suppression comes from the hierarchy of scales between the inverse 
^ " frequency scale and the Vainshtein radius, rather than the orbital radius of the pulsar 

system. Extensions of these results will have direct relevance to infrared modifications 
of gravity, such as massive gravity theories, which are known to exhibit a Vainshtein 
mechanism. Generalization to Galileon models with higher order interactions are likely 
to provide stronger constraints. 
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1 Introduction 

The discovery of cosmic acceleration has spurred a search for consistent modifications 
of gravity in the infrared, [1]. Theories in which the graviton acquires a mass, either 
softly as in the Dvali-Gabadadze-Porrati (DGP) model [2] or cascading gravity [3], or 
a hard mass as in the newly developed ghost-free models of massive gravity [4] or their 
bigravity [5] extensions, are a promising class of such modifications. These models 
commonly include light scalar degrees of freedom which arise from additional graviton 
helicity states. At the linear level, the scalar modes do not decouple in the massless 
limit, which is known as the van Dan, Veltman, Zakharov (vDVZ) discontinuity, [6]. 
Despite this, massive gravity models satisfy standard tests of gravity since the new 
scalars become strongly coupled near dense sources, which suppresses scalar gradients, 
and yields a force that is much smaller than the Newtonian one. This strong coupling 
effect is known as the Vainshtein mechanism, [7]. Far from compact sources, the scalars 
are weakly coupled and massive gravity theories make novel predictions that differ from 
Newtonian gravity in interesting ways. This suggests that IR modifications of gravity 
have extremely small effects on any but cosmological scales. 

Much of the phenomenology of the Vainshtein mechanism can be captured by 
considering the Galileon models [8]. Indeed the original Galileon model came from 
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considering the decoupling limit of the DGP model [9], and the structure of this scalar 
was generalized in [8] to include all interactions consistent with the Galilean symmetry 
7r — > n + c + v^x^, having a well-defined Cauchy problem. It was subsequently shown 
that the generic Galileon arises as the decoupling limit [10, 11] of the generic ghost- 
free massive gravity theory [4]. In an independent line of reasoning, the Galileon 
models have been viewed as scalar theories in their own right (independent of their 
graviton helicity-zero origin) and the most general covariant Galileon theories have 
been constructed [12-15]. In the decoupling limit in which the helicity-zero graviton 
mode is viewed as weak, all of these theories take on a similar structure and at least in 
some cases, the physics of the Vainshtein mechanism is qualitatively similar. For this 
reason, in this article we shall focus on the simplest cubic Galileon model, and will leave 
generalizations to subsequent work [16]. See Refs. [17-20] for studies investigating the 
Vainshtein mechanism directly in massive gravity and Refs. [21] for Fierz-Pauli massive 
gravity. For other potentially related observations of the Vainshtein mechanism, see 
Ref. [22-24]. 

In the simplest cubic Galileon (DGP) model, the predicted anomalous acceleration 
Aodgp is AanGp/aN ~ 10~ 15 , where a a? is the Newtonian acceleration, for a typical 
binary pulsar system of characteristic mass 3M and semi-major axis a ~ 10 -2 AU. 
We show that this ratio does not set the anomalous contribution to radiated power, 
which is a few orders of magnitude larger in known binary pulsar systems and could 
potentially be even stronger in slow, high eccentricity systems. The enhancement is 
due to the fact that the relevant length scale is set by the onset of the "far-field" 
region at radius fip 1 where f2p = 2n/Tp and Tp the orbital period (often denoted Pb 
in the literature.) For Tp = 8 hours, this is ~ 10 3 greater than the orbital radius. 
Furthermore, the system also radiates in the monopole and dipole channels since we 
are now radiating into a scalar. Monopole and dipole radiation exists because the 
scalar effectively violates the equivalence principle. These effects could potentially get 
more pronounced in future observations. 

Nevertheless, in the simplest model exhibiting the Vainshtein mechanism we con- 
sider in this paper, these enhancement are not sufficient to produce scalar gravitation- 
ally radiation effects which are within present observational limits from pulsar timing, 
and the Vainshtein mechanism, albeit slightly more subtle in this time-dependent sys- 
tem, is still very much alive. In a subsequent work we will show how higher order 
interactions in these different massive gravity/Galileon models are potentially much 
more strongly constrained by current observations [16]. 

The rest of this paper is organized as follows: In section 2 we review the cubic 
Galileon model and derive the general formalism to compute the radiated power in all 
generality in any multipole, taking care of the Vainshtein mechanism. The power is 
computed in two independent and equivalent ways using first an effective field theory 
approach and second a more conventional energy flux derivation. We then apply this 
formalism in the subsequent sections 3, 4 and 5 to compute the respective monopole, 
dipole and quadrupole radiation in the cubic Galileon model. Finally we compare these 
results with observations in section 6 and summarize our results. 
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2 Pulsar Radiation in Galileon Model 
2.1 Cubic Galileon 

Our starting point will be the so-called cubic Galileon model. This is the simplest 
model exhibiting the Vainshtein screening mechanism which was originally derived 
from the decoupling limit of the DGP model. In order to determine the emitted power 
it is sufficient to work in the so-called decoupling limit Mpi — > oo, in which the self- 
interactions of the helicity-two graviton are neglected. For this reason it will not be 
of importance whether the scalar field is viewed to arise from a Galileon model in 
its own right [8], DGP [2], massive gravity [4], or some other infrared modification. 
Nevertheless just to be explicit on our conventions, if we were dealing with the Galileon 
model the form of the action would be 



S 




2 4 V ; I 3A aV ' 



(2.1) 



where A is the strong coupling scale. In models of massive gravity, the scale A is 
associated to the graviton mass m by the relation A 3 = (m 2 Mp\). It is usual (although 
not essential depending on the context) to assume that the scale m is connected to the 
current Hubble scale m ~ H ~ 1.54 x 10 _33 eV, giving rise to a strong coupling scale 
A 3 ~ M Pl H 2 ~ (1000 km)" 3 . 

Writing the metric as g^ v = r]^ u + Mp^h^, we can work in the decoupling limit 
by sending Mp\ — > oo keeping A fixed. Assuming a conformal coupling to matter for 
7r, which is the case in all known theories of massive gravity (such as DGP, Cascading 
gravity or massive gravity), the action including matter is 

s = /d'x {-^(Shu - (i + ^) + ^hrr„ + JL «r 

(2.2) 

where (Sh)^ u = —}fAh ilv + . . . , and T is the trace of the stress-energy tensor. As we 
can see, the helicity-2 and -0 mode decouple in that case. The equations of motion for 
7T and h^y are then 

d, (-I^tt (1 + ^Dtt) + (dir) 2 ) = . 

Assuming a point source, T ^ = -M5^\x)5^5l, the background solut ion for 7r can be 
written simply if we introduce E(r) and take Vir(r) = fE(r), [25]. Then the equations 
of motion are solved by 



A 3 



9r 4 + — - 3r 2 

7T 



(2.4) 



where the Vainshtein radius r* associated with an object of mass M is 

M \ 1/3 If M N 1/3 



16m 2 J A V16Mpi 



(2.5) 
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For m ~ Hq, an object with mass 1M has r t = 3x 10 4 pc. We take the conventional 
branch solution E + which is trivial at infinity and bears no ghost-like instabilities. 



2.2 Perturbations 

The strong coupling effect is determined predominantly by the total mass of the system. 
We split the stress energy as T^ u = Tq V + ST 1 *" into a static 'background' part T^ v 
where we assume that the entire mass of the system is located at the center of mass 
(x = 0), and a perturbation 8T^ U which encodes the time dependent dynamics, which 
for slowly moving sources is 



5T? = - 



Mi5 {3 \x-Xi(t)) - M5 3 {x) 



i=l,2 



« , (2-6) 



where Mj is the mass of each companion and M = Mi + M2. 

Having decomposed the source in a background plus perturbations, we similarly 
split 7r = 7r + a/2/3 and express the quadratic Lagrangian for </>, 

£ -K i+ 4( £ ' + f))^4( i+ 4f)(^ 2 <"> 

2 V 3A' J V " v> V6Af P1 

This quadratic action will be sufficient to obtain our result for the power emission into 
scalar gravitational radiation. For pedagogical reasons we shall compute the power in 
two ways and demonstrate the consistency of the calculational results. 

2.3 Galileon Radiation 

2.3.1 Method I - Effective action 

In the effective action approach to the calculation of gravitational radiation [26], we 
first derive an effective action for the dynamics of the matter distribution by integrating 
out the graviton and scalar particles in the decoupling action (2.2), giving the non-local 
matter effective action expressed in terms of the Feynman propagator 

S cS = J d A xC M + J d 4 xd A x'5T(x)G F (x,x')5T(x') (2.8) 

+ usual helicity-two contributions from GR . 

Here we have used the fact that the field <fi can be expressed in terms of the Feynman 
propagator 

4>{x) = -^ — f d 4 x'G F {x,x')8T{x'), (2.9) 
V6Mpi J 

where we have defined the Feynman propagator via 

DG F (x,x') = i5\x -x') . (2.10) 
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The modified d'Alembertian is the one appropriate to the quadratic action (2.7) 

1 + 3^( S ' + f)) V ^' 
As usual, the Feynman propagator can be expressed in terms of Wightman functions 
G(x, x) = 6(t - t')W + (x, x') + 6{t' - t)W~(x, x) , (2.12) 

where 

poo 

W + {x,x') = J2 duut^rMKmir',^-^-^ , (2.13) 
em Jo 

and the ue m (r,Q)e~' luJt = ui{r)Yi m (VL)e~ luJt are a complete set of mode functions in 
spherical harmonic space which satisfy the homogenous equations of motion following 
from action (2.7). 

The time-averaged power emission P is 

p = -(^)= f (2-14) 



where f{oo) is determined from the imaginary part of the effective action 

2Im£ Gff f 00 

— = / dw/(w), (2.15) 

+ P JO 

and S e g is calculated over one period. We define the moments 

M lmn = -W P dtd 3 xu em (r, Q)e- innpt 5T . (2.16) 
i p Jo 

Taking the Fourier transform 

oo 

M lm = M lmn e mnpt , (2.17) 

n=— oo 

we have 

/( W )=™*T"E A f^'^{uj{t-t'))M, m {t)Mt m {t'). 
ZMpi+P^Jo J-oo 

From (2.14), the period-averaged power emission is then 

oo 

P = nfipl^&nnl 2 . (2-18) 

PI n=0 «m 
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Inside the strong coupling region r <C r*, the mode functions are approximately 

1/4 (y/% 



ue{r) = u 



J v * | — u 



(2.19) 



with 



(2£+l)/4for £>0 
-1/4 for £ = 



(2.20) 



and where we require the mode to be regular at r = 0. This boundary condition is 
actually ambiguous for the monopole, and the monopole solution is determined by 
requiring that the first derivative vanishes u' (r = 0) = so that this function is 
analytic in Cartesian coordinates. The constant u is fixed by imposing the correct 
normalization for the Feynman propagator 



DG F (x,x') = i5 4 (x-x') 



(2.21) 



Since = 0, we have 



nG F {x,x') 



2 /TV \ 3/2 



'2tx V r 

2 /TV \ 3/2 



6(t - t')d t W + {t, t') - 5{t' - t)d t W~{t, t') (2.22) 



2ix \r 



5(t-t') d t W + (t,t')-d t W-(t,t') 



We should therefore have 



\/f {j) 3/2 [W+foO - 9 t W-(t,t')] t=t , = -if 



(x — x') 



(2.23) 



(2.24) 



with 



2 /r*\ 
7r V r 



d t W + (t,t') 



t=f 



J-(-Y /2 J2 f°° duuu(r)u(r')Y etm (0, <t>)Y e %(8' , <f>') 



ir \ r 



Lm 







— IU 



- r 3/2 ^2, Y L™{ 



Lm 



duu (rr') 1 ^ 4 J v 



v * ■ — — ur J u + 



V3 
— our 
2 



-iu 



2 (rr' 



,1/4 



7T r 



3/2 3 



24«(r-r')^ v . 



7r 3 r 2 



Lm 



dqqJ u * (qr) J u ± (qr') 



-iu 2 r+\ —-S 3 (x — x') . 

V 7T 3 



(2.25) 



Consequently 



'2 /rv\ 
7r V r 



[d t W + (t,t') - d t W-(t,t% =tl = -2zu 2 n ^ 



5 3 (x-x'), (2.26) 
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and the Greens function has the correct normalization for 

Inside the strong coupling radius, the properly normalized mode functions are therefore 



We can use the WKB approximation to extend this solution out to r > r*, 

1/4 / x 

cos ( ^wr ) for u^ 1 <r<r 4 




1 1 



- cos (cor) for r>r. 

which corresponds to the correct field normalization well outside the Vainshtein regime, 
where the field is living in flat space-time. For rw < 1, we are well-inside the strong 
coupling region, and the mode functions can be expanded as 



, 9tt\ 1/4 1 frV/\ 3"/* 



128 / JK \ r* / 4 V 



r(l + i/*) 16r(2 + ^ 



(ru) 2 



.(2.30) 



Using these in (2.16) we can then derive the power emitted from (2.18), which we apply 
in various cases after giving an alternative and equivalent derivation of the radiated 
power. 

2.3.2 Method II - Energy flux 

Our second technique involves computing the energy flux for the ir field directly. In a 
diffeomorphism-invariant theory there is no local definition of the stress energy of the 
gravitational field. One can define a pseudo-tensor t^ u such that 

d^-g(T^ + ir)) = 0, (2.31) 

where T^ v is the stress energy of non-gravitational fields. The conservation condition 
is satisfied by t^ u = t^ L + t% u , where t^ L is the usual Landau-Lifshitz pseudotensor for 
the metric g^ u expanded to second order, and t% u is the stress energy of the n field 

t% = \ (d^ir - \g^A + ^^tkUDtt (2.32) 

(d^d<K) 2 d u -n + d v {d7i) 2 d^ - g.Md-n fd^) . 



4A 3 



Since the stress-energy t^ u splits into two decoupled terms, and the first contribution 
gives the standard Einstein gravity result, we need calculate only the additional flux 
from the 7r field. Expanding around the background solution 7r , the radial energy flux 
gives 

4 E" 



tl r = (l + ^g-ja^, (2.33) 
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up to irrelevant total derivative terms. The stress-energy is conserved, so we can find 
the total power emission by integrating the flux over a sphere of any convenient radius. 
The solution for an outgoing wave in the region r 3> Qp 1 is 

<f> = ^2(f>t mn (r)Y tm e inSll ' t , (2.34) 

Iran 

where 

4>£mn ~ 



y/Cs 0-lmn ^-jnQp f cj 1 dr (2 35) 



with the wave propagation speed c s = 2/\/3 in the strong coupling region. From this 
we infer the period-averaged power emission of 

oo oo 

P = J2 Yl n2 n 2 p\ai mn \ 2 = 2^^n 2 n 2 P \a, mn \ 2 . (2.36) 

n=— oo ira n=0 ira 

From the equation of motion for it is straightforward to derive the Wronskian identity 

Me mn 



lim 



AE 



3A 3 r 



y/EMpi 



from which we infer that 



|0£mn| 2 — n , 1 T7T I Me m n \ 2 , (2.37) 

6\n\\lp Mpj 

for n 7^ 0. On substitution into Eq. (2.36) this reproduces exactly the result of 
Eq. (2.18). Having derived the general formula for the power emitted by the Galileon in 
two independent ways, we now turn to computing the explicit low multipole radiation. 



3 Power Emission in Monopole 

To determine the power emission, we compute the multipole coefficients Mi mn and use 
Eq. (2.18). The Keplerian orbit of the two masses in the binary system is 

f(l-e 2 ) M 2 ,i 
ri - 2(t)= l + ecos^ ^' (3 ' 1} 

where e is the orbital eccentricity, and f the semi-major axis. In particular, near the 
pulsar, the zero mode function behaves as 

^-dW 1 -^)- (3 ' 2) 

where 

f3=[ — -0.69. 3.3 

' V 128 / r 3/4 
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3.1 Relativistic Corrections 

It transpires that as in GR, the leading order contribution to the monopole and dipole 
formula will vanish as a consequence of energy and momentum conservation despite 
the fact that the equation for the scalar fluctuations is strongly Lorentz violating in 
the Vainshtein regime. For this reason we need the mode functions to subleading order 
for the monopole and dipole to obtain a non-zero estimate for the power emitted in 
these modes. Given this fact we should work with the first relativistic corrections that 
arise to the trace of the stress energy. To first order in relativistic corrections we have 



5T 



=1,2 



1 2 

—vf 



$* + ...] 5 {3 \x - Xiit)) - M5 {3 \x) 



(3.4) 



where $j are the Newtonian potentials evaluated at the location of each source. How- 
ever as in GR, to second order in velocity, the conserved non-relativistic energy Enr is 
simply given as E N r = Yli(\MiV 2 ~ so that we can express M^vf + $;) = 

Mivf — E NR to that order in velocity. The velocity includes a radial and an- 
gular contribution, v 2 = f 2 (t) + (W^t)) 2 but using (3.1), we can check that the 
radial velocity is suppressed by two powers of the eccentricity compared to angu- 
lar velocity (our^i)) 2 . Since the eccentricity is of order e ~ 0.1 — 0.6 in realistic 
binary pulsar models, we can ignore their contribution, and use the approximation 
^ZiM^vf + $j) ~ £\ Mi(uri(t)) 2 — -Enr- To second order in the velocity, the trace 
of the stress-energy tensor thus include the relativistic corrections, 



5T = - 



^2 Mi (1 - K(t)) 2 ) 5®(x - S t (t)) — (M — E NR )6 {3 \x) 



i=l,2 



(3.5) 



Coupling the mode functions to this trace, therefore leads to the following contribution 
up to second order in our, 



d 3 xuo(r)ST 



P 



(our 3 ) 1 / 4 
5/3 



i=l,2 ^ 



4 



(cunit)) 2 ) — (M — Enr) 



(3.6) 



4(wr 



3)1/4 



Mi(ouri(t)) 2 + non-radiating contribution, (3.7) 



i=l,2 



where we ignore the piece going as Enr as it is conserved and thus cannot lead to any 
radiative contribution (time-independent contributions scale as 5 n o in the monopole 
moment M 00n and thus do not affect the power which goes as ^ n>0 n\M 00n \ 2 . This is 
to be expected since static sources cannot radiate). 



3.2 Power 

Defining the effective monopole mass 



M 



M 



M X M 2 + M 2 M 2 
W 2 



(3.8) 
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and using the expressions (3.6) and (3.1) in the relation (2.16) for the monopole mo- 
ment, with £ = m = gives 

Moon = J, ^LJ^ Mm (1 - e 2 ) 2 r ■ e ~ m "l . (3-9) 



4(a;r3)V4 Tp MV ' J (1 + e cos(fi P t)) 2 

The fact that it is the sub-leading term in the mode expansion that gives the main 
effect is a statement that the leading order contribution to the monopole radiation 
vanishes due to conservation of energy/mass, which is also why there is no monopole 
radiation in GR. The difference here is that because the source is ST and not the exactly 
conserved energy, there are relativistic corrections which lead to monopole radiation. 
Defining 



/?(e) = ^-(l-e 2 ) 2 / t— V2 , (3.10) 

Z7T J (1 + e cosx)^ 

the expression for the total power emitted via the monopole is hence 

-Pmonopolc = o^2 TiVlp \ M 00n \ 2 (3-H) 

P1 n>0 

K?S<S*f«w, (3, 2 , 

This integral /* f (e) is highly dominated by the low-n harmonics (till n = 10 as can be 
seen from Fig. 1). In principle one could include the sum of all the n's till a maximal 
n max = (fa;)" 1 at which point the approximation used in deriving the mode functions 
breaks down. However the contribution these extra harmonics is completely negligible. 

Performing this sum over |J* f (e)| 2 numerically for a typical eccentricity of 0.617 
corresponding to that of the Hulse- Taylor pulsar, we obtain, 

P , -74 ( " pf)4 Mll Cl 2 (3 13) 

^monopole -(A ^^3/2 P ' 1 J 

3.3 Comparison with GR 

The formula (3.12) is the analogue of the Peters-Mathews formula [27] for quadrupole 
gravitational radiation, 

32 G s Ml M 2 {Mi + M 2 )) (l + f e 2 + |e 4 ) 

-TPeters-Mathews _ 5 £ 5 f 5 f\_ g 2)7/2 ' 

which using Kepler's third law tip = G{M X + M 2 )/f 3 , c = 1, M|j = l/87rG can be 
re-expressed into current notation as 

A (l + Z3 e 2 , 37 4\ M 2 M 2 

P _ 4 V 1 1 24 e 1 96 C J /q s4 ^1^2 q2 (n 1 c\ 

Peters-Mathews - 5^ (1 _ ,2)7/2 l S W (Ml + M 2 ) 2 ' ( ^ 

Comparing our two results we learn one important thing: The suppression factor that 
arises due to the Vainshtein effect is given by 

Radiative Vainshtein suppression = 7— 77777- (3.16) 

(ilpr+y/ 1 
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This should be directly compared with the typical suppression from fifth forces 

/- \ 3/2 

Static (Fifth Force) Vainshtein suppression = I — ) . (3-17) 

\r*J 

In other words there is an enhancement of (fipr) 3//2 which explains why the 

Vainshtein effect, whilst still fully active, is nevertheless slightly less powerful in this 
fully time-dependent situation. 



4 Dipole Radiation 

We now turn to the dipole radiation for which the mode functions take the form in the 
small (ur) regime 

(cur' 



Ml(r) "(^( 7o ~ 7l(wr)2 )' ^ 



with 



7o= (^V /4 ^!_L^^o.4 and 7 i = 7o— r ( 7 / 4 ) ^o,04. ( 4 . 2 ) 
10 \128J 4 3 / 4 T(7/4) ' ro 16r(ll/4) v ; 

When computing the relativistic corrections to the trace of the energy momentum 
tensor, the situation is slightly different for the dipole than the monopole, and cannot 
easily be derived. Nevertheless, since the order of magnitude of the dipole is well below 
that of the monopole and the quadrupole, as it shall transpire, we can simply give an 
order of magnitude estimation and use, 

I d 3 x Ul (r)5T = -t-^ttj Yl M * M^W) - ci(^(t)) 3 ) , (4.3) 



where c and c\ are two constants with c\ ~ 70 + 71 ~ 0.4. 
4.1 Leading Contribution 

The dipole moment M^ n corresponding to the leading contribution in co(wrj(t)) is, 

where the angular position of the first object is taken to be localized at 6\ — tt/2, 4>\(t) = 
Qpt, while the second object is diametrically opposed 82 = tt/2, 02 (t) = Qpt + ir. Using 
the fact that Yi )0 (9 = 7r/2,0) = for any angle and 1x^1(^/2,0) = — lx,±i(0, + 71 ")? 
we therefore see, as in the case of the monopole that the leading contribution cancels, 
which this time is a consequence of angular momentum conservation. 

The leading order contribution to the dipole moment M^ n is the same as the 
dipole moment in GR, and this contribution vanishes for the same reason as in GR, as 
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angular momentum is conserved. We have already seen this with the monopole, where 
the leading contribution vanishes because of conservation of energy (mass in this non- 
relativistic limit). The difference between the monopole and the dipole is that in the 
present case the subleading contribution to the monopole is parametrically larger than 
anticipated, whereas for the dipole we shall see that it is much smaller. 

4.2 Subleading Corrections 

The next order moment M[H n in the expansion in the expansion in (cur) gives a non- 
trivial but small contribution 

(1) (erf (1-e 2 ) 3 /*> e-^ 

M lmn ^3)1/4 Tp J q (1 + £ cos(fipt ))3 

MtM^ ,tx _ , M 3 M 2 ,7r „ A , . 

_^y lm (-, n P t) + _L_ly lm (_ fip t + ,r) J , (4.5) 

which is only relevant for for m = ±1 and for systems with significantly different 
masses which is not the typical case of binary pulsar systems. Setting 

AM = 2 y 3 ^ , (4.6) 

we get 

The m = —1 channel is suppressed compared to the m — 1 as it picks up the opposite 
rotation direction. The main reason is that the m — 1 contribution goes e tm ^ and 
thus cancels the oscillating behaviour of e~ m< ^ for n — 1, whilst the m = — 1 never 
cancels the oscillation for positive n. An explicit computation, shows that P e= i m= _ 1 < 
10~ 2 P i= i m= i, and thus only focus on the m = 2 contribution in what follows. The 
power emitted through this channel is thus 

_c\ (gpff AM 2 - 



with 



r^ e ) = (1 - e 2 ) 3 ^— / — -dx . (4.9) 



PI n=0 



n 13/4 i-2tt e -i(n-l)x 



27r /n (1 + ecosx 



3 



Performing this sum over I n (e) numerically for a typical eccentricity of 0.617 corre- 
sponding to that of the Hulse- Taylor pulsar, we obtain, 



X>>.617)| 2 ^10 3 , (4.10) 



n=0 
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which gives an overall dipole power emitted of 



p -40 ^ 6 AM 'n 2 (A U ) 

Pdipoie ~ 40 (^rp ^ ^ I 4 - 11 ) 

which, for the Hulse- Taylor pulsar is roughly 9 orders of magnitude below the monopole 
radiation and is thus utterly negligible. Notice that the situation could potentially 
change ever so slightly for systems with a bigger mass difference between the compan- 
ions. 



5 Quadrupole Radiation 

We finish by estimating the quadrupole emitted by the Galileon. In that case the mode 
function is of the form 

"'M^S^ 1 -^ 2 )' (5 - 1) 

where 

3 9 / 8 7T 1 / 4 

2iV4r(9/4)- (5,2) 

In this case the leading contribution to the mode functions gives the dominant contri- 
bution to the power and we can thus ignore the subleading and relativistic corrections. 
The moment vanishes for m — ±1. Here again, out of the three other moments, 
m = —2, 0, +2, the leading contribution arises from the m = 2 one as it goes as e im ^ 
and thus cancels the oscillating behaviour of e~ m ^ for n = 2. For m = and —2 on the 
other hand this cancelation can never occur for n > and thus the contributions from 
these integrals is always suppressed compared to the m = 2. An explicit computation, 
shows that P e=2)m=0 < lQ- 2 P i=2m=2 and P t=2 ,m=-2 < 10~ 3 Pi=2,m=2, and thus only 
focus on the m = 2 contribution in what follows. In that case 

M, M ,M = -i^^g^tf (e) , (5.3) 



with 



and 



M * = M3/2 ' > ( 5 - 4 ) 



7/4 f 2, 



/«( e ) = (l_e 2 ) 3 / 2 ^— / — -dx. (5.5) 

nK J v ; 2tt J (1 + ecosx) 3 / 2 v ; 



The power emitted through this channel is thus 

r2 00 



n=0 
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Figure 1. Numerical evaluation of the monopole. dipole and quadrupole integrals iff (e), 
In(e) and I„ (e) expressed in (3.10), (4.9) and (5.5) respectively for the eccentricity of the 
Hulse- Taylor pulsar 1913+16, e = 0.617. The right hand side represents the numerical 
evaluation of the sums Y^N=o\^n i e )\ 2 ■> Sw=o \^n( e )\ 2 an d Yln=o\^( e )\ 2 ^ w hich shows a 
rapid convergence as N — > oo. 

Once again, the sum is dominated by the low-frequency harmonics n < 15 but the sum 
can be performed numerically for up arbitrarily large values, (see Fig. 1). Fortunately 
it converges rapidly, and for the Hulse- Taylor pulsar, we obtain 

oo 

^|/«(0.617)| 2 ^18, (5.7) 

n=0 

which leads to a power emitted roughly 10 times larger than the monopole. 
6 Discussion 

Pulsars in DNS binaries are used to measure gravitational-radiation effects. Timing 
measurements of these systems are free from contamination due to tidal effects or 
accretion from a stellar companion. The archetypal such system is the Taylor-Hulse 
pulsar 1913+16. We have summarized the orbital parameters of four known DNS 
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pulsars (A to D) and one pulsar- white dwarf binary (E) in Table 1, [28], all of which 
have measured orbital period derivatives Tp, which agree with GR. 

The orbital period derivatives Tp are given in terms of the non-relativistic energy 
i?NR and the power emitted by the relation 

t P = ^T P 1 d -^ NR (g ]_) 

3 -Enr dt 

where the system's non-relativistic energy is 

1 M X M 2 ( Qp \ 1/3 

NR = 0/0-^/3 JW WW ' I 6 ' 2 ) 



2(8tt) 2 / 3 Mpi V mm pi 

The monopole and quadrupole Galileon radiation are close, but in the explicit exam- 
ples presented here, the Galileon quadrupole always give the largest contribution. In 
summary, in the simplest model, the Galileon radiation is at least 7 orders of magni- 
tude below that of GR (7 orders of magnitude for pulsar C and E and 8 for the other 
DNS pulsars), which 6 orders of magnitude below the current precision, [36], when 
considering a parameter m = 1.54 x 10 _33 eV, or equivalently for a strong coupling 
scale A ~ 10~ 13 eV ~ ( 1000km) . The best precision is still the Hulse- Taylor pulsar 
and pulsars C and E. Even though the double pulsar D has a good precision a, its low 
eccentricity make it not the best candidate to probe the Vainshtein mechanism. 

As can be seen from the expression (6.1) for the orbital period derivative, using 
the power emitted in the monopole, (3.12), the dipole (4.8) or the quadrupole, (5.6), 
tp scales directly as r* 3 ^ 2 ~ m ~ A 3 / 2 . So in order for the scalar field 7r to have an 
effect at all on current or upcoming binary pulsar timing observations, the parameter 
m should be enhanced by roughly 6 orders of magnitudes. Binary Pulsar timing thus 
put a rough bound of m < 10~ 27 eV or A < 10~ 9 eV. Compared to solar system tests [? 
], these bounds are not competitive for the cubic Galileon interactions, but could be 
enhanced when considering higher interactions, [16]. 

Present constraints are limited by the sample of DNS pulsar systems. Ideal DNS 
binary pulsars have long periods, high eccentricities, and are located nearby so that 
the kinematic corrections are reduced. Over time, measurements of the orbital period 
derivative become more precise as more data is gathered. The ultimate precision is 
limited by the uncertainty in the relative acceleration between the sun and the pulsar 
system, which must be included to obtain the intrinsic orbital period derivative from 
the apparent one. 

We should note that whilst previous authors have considered constraints on the 
mass of the graviton in binary pulsars, in particular see [38], these authors do not take 
account of the Vainshtein mechanism, working only in the linearized Fierz-Pauli theory, 
and they further utilize an incorrect expression for the stress energy of gravitational 
radiation which accounts only for the helicity two component. As such the constraints 
obtained there are not appropriate for consistent Lorentz invariant theories of massive 
gravity as considered in [4]. They may however be relevant to theories in which there 
is no propagating helicity zero mode. 
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In this paper, we have shown the successful implementation of the Vainshtein 
mechanism in a fully time- dependent setup. This explains how a conformally coupled 
scalar field, or the helicity-0 mode of a massive graviton can evade typical tests of GR 
and particularly the well constrained orbital period decay in binary pulsar systems. 
Nevertheless, despite the presence of an active Vainshtein mechanism, we show that 
the suppression is less than naively anticipated from purely static systems for several 
reasons: 

• The suppression factor that arises due to Vainshtein effect is going like (fipr*) -3 / 2 
rather than (f/r*) 3 / 2 as is the case in static spherically symmetric configurations. 

• Furthermore, even though the leading contributions from the monopole and 
dipole radiation vanish from energy and angular momentum conservation, the 
sub-leading (relativistic) contributions are non-negligible, and in the case of the 
monopole can be comparable to the Galileon quadrupole radiation. 

• Finally, we have focused here on the simplest realization of the Vainshtein mech- 
anism, namely within the context of the cubic Galileon model. However higher 
order Galileon interactions could potentially lead to an additional enhancement 
of the Galileon radiation and ought to be studied in their own right, [16]. 
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